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Abstract 

Elastic quantum bound-state reflection from a hard-wall boundary provides direct information 
regarding the structure and compressibility of quantum bound states. We discuss elastic quantum 
£T) bound-state reflection and derive a general theory for elastic reflection of shallow dimers from 

hard-wall surfaces using effective field theory. We show that there is a small expansion parameter 

<N 

for analytic calculations of the reflection scattering length. We present a calculation up to second 

a 

order in the effective Hamiltonian in one, two, and three dimensions. We also provide numerical 
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lattice results for all three cases as a comparison with our effective field theory results. Finally, 
we provide an analysis of the compressibility of the alpha particle confined to a cubic lattice with 
vanishing Dirichlet boundaries. 
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I. INTRODUCTION 

The study of confined quantum bound states has interesting applications in many areas 
including nuclear structure calculations, experimental ultracold atomic physics, and quan- 
tum dots and wells. Elastic reflection off of a hard wall allows one to gain insight into the 
structure of a bound state as well as how it behaves when confined or compressed. In the 
case of quantum dots and wells, the bound state is comprised of an electron and hole con- 
fined inside a two- or three- dimensional nanoscale semiconductor structure. Varying the 
geometry of the structure allows increased control of current tunneling as well as photon 
absorption and emission. See for example Ref. [IrH]. The theory developed in this paper 
gives a universal result for the energy of Wannier excitons in direct-band gap semiconductors 
as a function of the binding energy, effective masses, and the geometry of the nanostructure. 

In the case of ultra-cold atomic experiments, a quantum well is produced by using lasers 
to create repulsive boundaries, thereby confining the atoms on an optical lattice. Similar 
ideas have been proposed for quantum billiards systems j5]. Since they add no dimensionful 
scale to the problem, hard-wall boundaries provide a probe of universal physics at large scat- 
tering lengths. For example the scattering length for particle-particle scattering is directly 
proportional to the scattering length for dimer-wall reflection. In this paper we calculate 
this universal proportionality constant in one, two and three spatial dimensions. 

It is not possible, experimentally, to construct a hard-wall boundary for protons and 
neutrons, so it may seem that the current discussion has no direct connection with nuclear 
physics. A similar critique could be made of Liischer's analysis of periodic boundaries in finite 
cubic volumes [5J [7] . However, Liischer's analysis now provides the theoretical framework for 
numerous calculations in lattice quantum chromo dynamics [8HTU] . While it is not possible to 
construct a hard-wall boundary for protons and neutrons in nuclear physics experiments, it 
is possible to do so in ab initio numerical calculations. In some cases such calculations can be 
used to verify phenomena observed in experiments, and in other cases the calculations probe 
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new physics inaccessible in the laboratory. As with temperature and chemical potential, 
boundaries can provide tunable control parameters for use in such calculations. Hard-wall 
boundaries can and have been incorporated into ab initio numerical calculations [TTJ [T2] . 
Hard-wall boundaries are currently being used to perform nuclear lattice calculations that 
probe properties such as structure and elastic deformation of nuclei. The method is a useful 
compliment to Luscher's periodic boundary analysis for bound states and scattering states. 
By varying the boundary conditions one can calculate quantities such as nuclear radii and 
quantities analogous to bulk and shear modulus, which can then be compared with bulk 
modulus estimates inferred from observed monopole resonance energies. 

Recently there has been a great deal of interest in alpha-particle clusters confined inside 
nuclei such as carbon-12 [T3HT5] . Very recently ab initio lattice effective field theory calcu- 
lations of carbon-12 have given the energies of the ground state, the excited spin-2 state, 
and, for the first time, the Hoyle state responsible for the formation of carbon in stellar 
environments [IB]. Furthermore, these calculations indicate the presence of compressed, cor- 
related alpha clusters. Since there are no low-energy resonances of the alpha particle and 
little experimental data on the compression of alpha particles, a better understanding of 
this phenomenon would prove valuable. In this paper we analyze the compressibility of the 
alpha particle. 

In this paper we discuss the elastic reflection of quantum bound states off of hard-wall 
boundaries. We begin with a short presentation of the general case, where the number of 
dimensions and number of constituent particles is completely arbitrary. We then develop 
a general theory from the principles of effective field theory for the specific case of shallow 
two-body quantum bound state reflection in one, two, and three dimensions. Our main 
result is a derivation of the phase shift due to the scattering of a two-particle bound state 
on a hard wall in the adiabatic limit and up to second order in an expansion of the effective 
Hamiltonian. The effective field theory results are presented for the cases of one, two, and 
three dimensions for arbitrary mass ratio 777,2/7711 and compared to numerical results. Finally, 
we provide an analysis of the compressibility of the alpha particle, summarize, and consider 
further work to be done. A partial summary of our results was presented in a previous 
letter publication [17] . Here we present the full details of our calculations. We note that 
very recently there has been work on two body systems in a finite volume with periodic 
boundaries [T5J [TJJ] as well as two body systems with Neumann boundary conditions [2D] . 
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FIG. 1: Sketch showing the distance X from the wall to the center of mass of the bound state and 
the standing wave solution 

II. FORMALISM 

We consider a non-relativistic bound state in d dimensions, with total mass M, and 
arbitrary number of constituent particles. The bound state scatters elastically off of a hard- 
wall boundary, implemented as a vanishing Dirichlet boundary condition. We let X be the 
distance from the wall to the center of mass of the bound state. We work in the inertial 
frame where the only non-zero component of the center-of-mass momentum is the component 
perpendicular to the hard-wall boundary. Then we construct a standing wave solution with 
center-of-mass momenta ±p perpendicular to the wall with reflection phase shift, 5{p), and 
reflection radius, R(p) = —S(p)/p. The reflection radius is a measure of the distance between 
the wall and the closest node of the asymptotic standing wave, ty p (X) oc sin [pX + 6(p)]. 
This is shown in Fig. [T] 

Reflection in one dimension is analogous to S'-wave scattering in three dimensions and so 
we use the effective range expansion, 

p cot 8{p) = — — + \r R p 2 - V R p 4 H , (1) 

Or Z 

where clr is the scattering length, is the effective range, and Vr is the shape parameter. 
Notice that at threshold, clr is equal to the reflection radius, clr = \im p ^o R(p). For a 
completely rigid bound state, R(p) = clr for all p. However, we expect that the bound state 
will not be perfectly rigid and so it will compress more when the collision energy increases. 
This means that the reflection radius, R(p), will decrease with increasing bound state center- 
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of-mass momentum, p. The rate of decrease measures the compressibility of the bound state 
under unilateral stress. 

We now consider the bound state confined to a <i-dimensional cubic box with length L 
and hard wall boundaries on all sides. In the limit of large L, we can approximate the 
wavefunction as a product of standing waves along each coordinate axis. Each standing 
wave has half-wavelength equal to L — 2R [p(L)], and so the reflection radius can be related 
to the ground state confinement energy of the bound state as follows: 

E(L) = ^§^ll + 0(L^)], (2) 

where 

P{L) = L-2 % R[p(L)Y (3) 
This relation can be used to determine the reflection radius as a function of the center- 
of-mass momentum. The O (L~ 2 ) relative error in Eq.Q arises because the different coor- 
dinates for d > 1 cannot be exactly separated as a product of standing waves. This is due 
to the effects of double-wall collision near the wall intersections. By wall intersections we 
mean the corners of a two-dimensional square or the edges of a three-dimensional cube. In 
each case, the codimension of the wall intersections is two, and this explains the O {L~ 2 ) 
relative error. In one dimension, however, there are no wall intersections, and so the error 
is exponentially small in L. 



A. General Two-body Bound States 

We now consider the more specific case of a bound state with zero orbital angular mo- 
mentum consisting of two distinguishable particles with masses mi and m 2 . The reduced 
mass, /i, is defined in the usual way, ji = (m^ 1 + We let Er be the infinite volume 

binding energy of the dimer, and ^ be the binding momentum defined by the relation 
Eb = — K fi/(2/i). dB is the 5-wave scattering length for shallow-binding particle-particle 
scattering. In the shallow-binding limit, where as is much larger than the range of the 
interaction Kb = g^ 1 , we can neglect the short-distance physics. So, the reflection phase 
shift is a universal function of the dimensionless ratio p/k,b- 111 this paper we present the 
form of this universal function. 
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B. Shallow Dimer Reflection 



Next we use the principles of effective field theory to calculate the reflection scattering 
length for a shallow dimer. While effective field theory is a well established method, in the 
current case we are dealing with a non-homogenous system for which there is no immedi- 
ately obvious small expansion parameter. In the soft scattering limit, we use an adiabatic 
approximation for the center-of-mass motion of the system. We then perform an asymptotic 
expansion to account for the long-distance physics. The result is an expansion for KbCLr 
with convergence controlled by an expansion parameter e~ 2KBaR . 

We let r\ and r 2 be the coordinates of the two constituent particles and assume an 
attractive short-range interaction given by the Hamiltonian, 



H —^-^ +c ^-^' (4) 



where 5^ is a regulated d- dimensional delta function. The coefficient Cb is tuned to produce 
a bound state with energy Eb at infinite volume. We take r = r\ — r*2 to be the relative 
separation between the particles. For any fixed center-of-mass coordinate, the part of the 
Hamiltonian that is only dependent on the relative coordinate, r, is given by 

H Tel = -^V 2 r + C B S^(^. (5) 

1. The Effective Hamiltonian 

Let Ek be the kinetic energy of the moving dimer. To calculate the reflection scattering 
length it suffices to consider dimer- wall scattering in the limit Ex <C \Eb\- In this low-energy 
limit we perform an adiabatic expansion for the center-of-mass motion. This technique is 
conceptually similar to the adiabatic hyperspherical approximation [211- 23j. 

We let X be the distance from the wall to the center of mass of the dimer. We label 
the coordinate axes 1, 2,...,d, and take the d th dimension to be perpendicular to the wall. 
For fixed X the hard- wall boundary at {r\)d = gives a minimum value for the relative 
coordinate r d , 

M 

rr - —X. (6) 
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We now define a;_ as twice this distance, 

,_<*) = A. (7) 
m 2 

So X- (X) is twice the distance from the wall to r 2 when fi touches the wall. In other words 
^x(f) must vanish when equals —x_(X)/2. The factor of two simplifies the expansion to 
be derived later. Similarly the hard-wall boundary for r 2 gives a maximum value for r d , 

M 

rT x = —X. (8) 
mi 



As in the previous case, we define 



X+ (X) = —X. (9) 
mi 



So x + (X) is twice the distance from the wall to fi when ri touches the wall. In other words 
ipxif) must vanish when r d equals x + (X)/2. Once again the factor of two has been added 
to simplify the expansion to be derived later. 
We now define the ci-dimensional vectors 

f_(X) = (0,---,0,x_(X)), (10) 

f + (X) = (0,.. - ,0,x + (X)). (11) 
The magnitudes of the vectors are 

r_(X) = |f_(X)| = -x_pO, (12) 

r + (X) = \f + (X)\=x + (X). (13) 

Let E K be the kinetic energy of the moving dimer. For each X it is only necessary to keep 
the ground state of if re i satisfying the hard-wall boundary condition. This is due to the fact 
that the contribution of the excited states are suppressed by powers of Ek/ \Eb\ and hence 
cannot contribute to the reflection scattering length, but only to higher-order coefficients in 
the effective range expansion. 

For each X we take tyx{r) to be the normalized ground state wavefunction of H TC \ sat- 
isfying the hard- wall boundary constraint. Then the eigenstates of the full Hamiltonian 
are 



MX))®\i; x (r))+0(PyKl), 
8 



(14) 



where P is the center-of-mass momentum of the dimer in the direction perpendicular to the 
wall. The low-energy effective Hamiltonian is 

H "--wi& + v{x)+T{x) - (15) 

where V(X) is the adiabatic potential, 

V{X) = (il> x \H td \il> x ), (16) 

and T(X) is the diagonal adiabatic correction, T(X) = — ^ (ip x \ ^2 \i>x) ■ 

Notice that the cross-term resulting from one derivative with respect to X vanishes be- 
cause of the fixed normalization of ip x . We use lattice regularization to deal with the 
continuum limit singularity in the delta function for d > 1. 

We define the local X-dependent energy, E^_^ + , and binding momentum, «f_ i f+,, m terms 
of the adiabatic potential 

k 2 

V(X) = 0j x \ H rcl |Vx> = E f+if _ = — y^, (17) 

for the vanishing Dirichlet boundary condition. For large X we can generate an asymptotic 
expansion for the local X-dependent binding momentum, K?_f + , V(X), and T(X) in powers 
of e~ Ksr +( x ) and e~ KBr ^ ( - x \ The boundary conditions are enforced by applying the method 
of images. 

We consider two infinite periodic chains of delta functions. The first chain includes 
delta functions centered at multiples of r + — r_. The second chain includes a delta function 
centered at r + plus multiples of r + — r_ . We now construct a wavefunction which is the lowest 
energy eigenstate for this potential with alternating signs at the centers of the neighboring 
delta functions. For a picture of the wavefunction see Fig. |2j 

For this wavefunction we use an ansatz that is a superposition of d— dimensional Yukawa 
functions. 

Ipr- ,r+ (r) = A?_ f + (fir?_ f + (f ) , (18) 

where 

<f>?-,?+{r) = yd,o(«f_,f +) |^|) - yd,o |r - f_|) - y d) o(Kfr_ >P+ , \f-f+\) 

+ yd,o(Kr-,? + , \r — r + +f-\) + yd,o(^r_,r + , r + r + — r_) H , (19) 
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physical region 




-2x + +2x. 



FIG. 2: Sketch of the wavefunction derived from the method of images 



Here 



yd,o{K,r) 



dp e 



-ip-r 



(2ix) d P 2 + k 2 ' 

is the generalized d dimensional Yukawa function, and we also define 



(20) 



r) = d d r" yd,o{^, \r'\)y d ,o(K n , \r' - f|), 



(21) 



as the full overlap integral of two Yukawa functions in d dimensions whose centers are 
separated by a spatial distance r. In Eq.(18) A?_^ + is a function that normalizes Vv_ ,r+ (^0 > 
and we find that 



d d p 

E 



1 



-E 



EL _ E~ - 2- 
d d k e~ i P'^~^ +1 ^ r++ ^ r -^ 



B 



(2tt) c 



2- - E~ - 

2fi ^r-,r+ 



-E 



d d k e - i P-i-ir++jr-) 



(2*) 



2ft 



Ef_ r 4 



(22) 



where j is summed over integer values. 

At first order in powers of e~ KBr +( x ^ and e -Ksr-< - x - ) we get a correction to the binding 
momentum 



K 



r_,r+(l) 



d d p 

(27T) d 



a 

8kb 



I 



d d p 1 
(2n) d P 2 +k|_ 



(23) 
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The first-order correction to the adiabatic potential is 



-Vdfi («B,r + ) 



and the first order diagonal adiabatic correction is 

T W(X) = M 

2mlf d (K B ,K B ,0)' 

2mlf d (K B ,K B ,0) yd '° KB,r ~ 
So the first order correction to the effective Hamiltonian is 



V\X) + T X (X) = + 

M 

+ 



M 



2mf/ d (« B ,K B ,0) 
2/d,o («s,r_) . 



Vd,o{K B ,r + ) 



2mlf d (K B ,K B ,0)' 
At second order the correction to the local binding momentum is 



r_,r+(2) 



afel/d.o (kb, 0) 



-2yd,o (^b, r+ + r_ 



yd,o(«B,7- + )+yd, (KB,r_) g 

^y d ,o {k b , o) 9k b 



[2/d,o(«B,f*+) +y d fl(K B ,1 



[y d fi(K B ,r + ) +y dfi (K B ,r_)} 2 d 2 



8 
8 k b 



Vdfl («B,0) 



9«| 



2/d,o («b,0) 



The correction to the adiabatic potential is 



The diagonal adiabatic correction contains a number of terms which we unite as 

T ( 2)(x) = T m {x) + T (2) (x) + T <g) (x) 



The first term in Eq.(29) is 



T( 2 Vm - fd( K B,K B ,0) f dA?_ M) 
Ha)W~ 2M \ dX 

where the ellipses indicate terms at third order or higher. Here 



+ 
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+2/ d r+ + r_) H ] 1/2 , (30) 

arises from the normalization of the relative coordinate wavefunction. Once again the ellipses 
indicate third order or higher terms. For more details see Appendix |Aj The derivative of 
the normalization function, at first order, with respect to X is 



if*_ ,r+ (1) 

dX 

2 [/ d (^s,0)] 3/2 
The derivative of the first order correction to the binding energy with respect to X is 

<9/V_,r+(l) 



(31) 



dX 



afcfd,o(«B,0) 



(32) 



The second term in Eq.(29), due to two derivatives with respect to the spatial coordinates 

ir_,r+(i) and /v_,r+(i), is 



without any contribution from the first-order corrections A f 



M 



+ 



HJd{KB,KB,V) 



(33) 



The third term in Eq.(29), proportional to «fi,f+(i) and is 



^r_,r+(l) 



3 



2mf / d (k, «, 0) 9k 



+ 



r_,r+(l) 



<9 



2m|/«j (k, k, 0) 8k 



[Vd,0 («,7*_)] 



+ 



.4 



r_ ,r+ (1) 



[/d(«B,«B, 



,1/2 



.1/ . M . 
— 2 3/d,o r +) H 2^,0 \^Bi r_ 



k=kb 



m 



(34) 



where y^i («, r) , (k, r) , and y<2,3 ( K , r ) are the generalized Yukawa functions defined later, 



in Eqs. (43), (52), and (59), for d = 1, 2, and 3. For the terms in the second line of Eq.(29) 
we find 



k b 



8k. 



f— ,r+ (1) 



(Dl) 



"T.i/ d (K B ,K B ,0) dX 

+y?/d,2 (kb, r+) + ^,3 



-3/d.i l«B,r +> 



(35) 
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m 2 fd(^B,^B,0) dX [8 
r _ 3 1 

+yyd, 2 («b, r_) + -y d , 3 («b, r_) L (36) 



9«r_,r + (l) X 2 



T (2) / Y N _ V dX ) f o d Vdfi ( n x 

+3/t ^!f ( ^' 0) - K2 Sf (kb ' 0) } ' (39) 



and the effective Hamiltonian is 



^ = - " j«-'*> ~ k'-^ + Ti2,(x) - (40) 

where /v_,r + (i), «f_,f+(2), and T^ 2 \X) are as defined above. 
III. RESULTS 

In this section we present the results of our effective field theory expansions in one, two, 
and three dimensions. Using the definitions given below for the generalized Yukawa function, 
full two Yukawa function overlap integral, and effective Hamiltonia terms from the previous 
section, we present the expressions for the effective Hamiltonian in one, two, and three 
dimensions. In all cases at zeroth order we start with the infinite volume result, 

#(°) - 1 d l_ + V (°) + T (0) (41) 

where 

-(0) , rp(0) 



V£>+T$> = E B . (42) 
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A. Results in One Dimension 

The generalized Yukawa function in one dimension is 

e~ KT 

2/i,o = — , 

and the full two Yukawa function overlap integral in one dimension is 



e 

h («, K,r) = — (1 + rer) . 



L Firsi Order Effective Hamiltonian in One Dimension 

At first order the effective Hamiltonian for the one-dimensional case is 



ttW _ et(0) , T/ (i) , T (i) 



where the effective potential is 



y(i) (x) + T (D (x) = 4^! 

mim 2 



e ~K B r + (X) e -K B r-(X) 
+ 

mi m,2 



2. Second Order Effective Hamiltonian in One Dimension 

At second order in one dimension the effective Hamiltonian is 



H$ = H V + T W(X) + VW(X), 



where 



and 



v(2){x) = ^ii ( vi2)(x) + ^ 2)(x) + ^ 2 ) )(x) ) ■ 

vi 2 \x)=e- 2K ^(2K B r ± -l), 
V$(X) = e~ KB ^ +r -^ (n B r + + k,bT— - 3) 
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T^\X) = T$JX) + T{l\(X) + 2%>(X) 




-i.o L 

FIG. 3: Plot of the effective potential at zeroth, first, and second order in one spatial dimension 
for mass ratio m-zjm\ = 4. 

Fig. [3] shows a plot of the effective potential, at zeroth, first, and second order, in 
one spatial dimension for mass ratio m^/wi = 4. The effective potential is plotted in 
dimensionless units of Xkb along the horizontal axis and M l ff along the vertical axis. 

Kg 

We see that the effective potential crosses zero at Xkb ~ 0.5, and the convergence for the 
expansion is good for Xkb ^ 0.5. This part of the potential gives the dominant contribution 
to the reflection scattering length. 

B. Results in Two Dimensions 

The generalized Yukawa function in two dimensions is 

2/2,0 («,?") = ^-Kq(kt), (52) 

Z7T 
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and the full two Yukawa function overlap integral in two dimensions is 



f 2 (k, k, r) = - — K^nr), (53) 



where K n is the modified Bessel function of the 2 kind. 

1. First Order Effective Hamiltonian in Two Dimensions 
In two dimensions the first order effective Hamiltonian is 



H$ = H® + (± + *lK [ KB r + (X)} 



+ [- + — 2 )k b K [k b t.{X)} ) (54) 



where the adiabatic potential is 



ji m\ 



7(1) (X) = [K (K B r + ) + K (n B r^] , (55) 



and the diagonal adiabatic correction is 



T«(X) = ^K Q (K B r + ) + £M Ko (K B r-). (56) 
mi mi 



2. Second Order Effective Hamiltonian in Two Dimensions 

At second order in two dimensions we find the adiabatic potential 



The diagonal adiabatic correction is 



t^\x) = t$(x) + rgpo + jM(x) 

+ T^X) + T$ 2) (X) + T^{X) + T^(X) + T$(X). (58) 
For the definitions of /V-.r+ti), K r_,r+(2), Vjf)> ^(B)> ^(C)> -^(di-D4)' an< ^ see Appendix 

o 
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FIG. 4: Plot of the effective potential at zeroth, first, and second order in two spatial dimensions 
for mass ratio mijmx = 4. 



m 2 /mi = 4. The effective potential is plotted in dimensionless units of Xkb along the 
horizontal axis and ^£ii along the vertical axis. Again the effective potential crosses zero 
at Xkb ~ 0.5, and the convergence of the expansion is good for larger values of Xkb- 

C. Results in Three Dimensions 

The generalized Yukawa function in three dimensions is 



Fig. [4] shows a plot of the effective potential in two spatial dimensions for mass ratio 




— KT 



(59) 



and the full two Yukawa function overlap integral in three dimensions is 




(60) 



1. First Order Effective Hamiltonian in Three Dimensions 



In three dimensions at first order the effective Hamiltonian is 
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H$ = H^ + TV(X) + V^(X), (61) 
where the effective potential is 

T (D(X) + V (1 \X) = KbM v \ e -^ + {x) + e - KB r_ ( x)i (62) 



2. Second Order Effective Hamiltonian in Three Dimensions 



At second order in three dimensions we find that the adiabatic potential is 

V {2 \X) = -J^ +(2) - ^4_A(D- ( 63 ) 
For the diagonal adiabatic correction, we find 



T^(X) = T{%{X) + T™(X) + 2${X) 

+ ^{X) + T$ 2) (X) + T^X) + T^(X) + T$(X), (64) 

For the definitions of Kf_,r + (i)> K r-,r + (2), ^m)> -^(bv ^(C)> ^(m-D4)> an< ^ "^(s)' see Appendix 

El 

Fig. [5] shows a plot of the effective potential in three spatial dimensions for mass ratio 
7712/777.1 = 4. As in the one and two dimensional cases, the effective potential is plotted in 
dimensionless units of Xkb along the horizontal axis and ^ Ve 2 S! along the vertical axis. The 

Kg 

effective potential crosses the horizontal axis at Xkb ~ 0.4, and the convergence of the 
expansion is good for Xk b > 0.4. 



IV. SCATTERING LENGTH AND REFLECTION RADIUS 

From these effective potentials it is straightforward to compute the reflection scattering 
length up to second order. This process can be carried forward to any order. The net 
result is an expansion with an expansion parameter of size e~ KB ' r± ^ aR > < e~ 2KBaR . The larger 
kbclr, the faster the convergence of the expansion. First- and second-order results for the 
one-dimensional system are shown in Fig. [6j Results for the two-dimensional system are 
shown in Fig. [7j and results for the three-dimensional system are shown in Fig. |8j In each 
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-1.0 



FIG. 5: Plot of the effective potential at zeroth, first, and second order in three spatial dimensions 
for mass ratio 1712/ mi = 4. 

case the agreement with lattice results is consistent with third-order corrections of size < 
A. Numerical Results 

We have calculated the dimer-wall reflection phase shift using Hamiltonian lattice meth- 
ods in one, two, and three dimensions. We implement the two-particle interaction as an 
attractive point-like interaction. The lattice Hamiltonian in d dimensions is given by 



H = — J2 Yl a \ (™) ai 

'=1,2 n 

d 

\a\ (ft) a, (n + Ij + a\ (ft) (n — I 



1 

2m 



i=l,2 n «=1 



ft) a\ (ft) a} 2 (n) a,2 (ft) , 



(65) 



where a\ and at are the creation and annihilation operators for particle i, n labels the lattice 
site, and I is the unit vector in the I th direction. We take parameters and operators to be 
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dimensionless by multiplying physical values by the appropriate power of the lattice spacing. 

Two parallel hard- wall boundaries are spaced a distance L z apart. We impose periodic 
boundary conditions in the perpendicular directions with length L, and consider dimer states 
where the momentum is perpendicular to the wall. We calculate the confinement energy as a 
function of L z using sparse-matrix eigenvector methods and from this calculation determine 
the reflection phase shift using Eq. g). For various mass ratios of the two particles, we repeat 
our calculations using successively smaller lattice spacings. This allows us to extrapolate the 
reflection radius and dimer kinetic energy to the continuum limit and determine universal 
results in the shallow-binding limit. For the two- and three-dimensional cases we also perform 
infinite volume extrapolations in the dimension perpendicular to the wall. 



The coefficient of the delta function interaction in Eq. (65) is c. In order to take the 
continuum limit in one dimension, we fix L z \c\ and extrapolate both Rkb and Ek/\Eb\ for 
three different values of L. The values of L z \c\ we consider are L z \c\ = 8, 10, 20. For L z \c\ = 8, 
we calculate the ground state values for Rk b and E K /\E B \ for L z = 80,160,240. For 
L z \c\ = 10, we consider the ground state and first two excited states for L z = 100,200,300. 
For L z \c\ = 20, we consider the ground state for L z = 200,300,400. 

In two and three dimensions, we first take the infinite volume limit in the perpendicular 
directions by extrapolating L — > oo. We then extrapolate to the continuum limit using 
several values for L z while keeping k b L z fixed, where Kb is the dimer binding momentum. 
In two dimensions we consider the ground state and first excited state for L z k,b = 6.8, 10.2. 
For L z k b = 6.8 we perform calculations using L z = 40, 60, 80, 100, and for L z k b = 10.2 we 
consider L z = 60, 90, 120, 150. 

In three dimensions we consider the ground state and first excited state for L z k,b = 
4.9,6.0. For L z kb = 4.9 we perform calculations using L z = 20,30,40, and for L z kb = 6.0 
we consider L z = 24,36,48. For a summary of the continuum extrapolations performed, 
including the lattice volumes used, see Appendix |E} 



1. Numerical Results in One Dimension 

Fig. [6] shows the one-dimensional results in the shallow-binding limit, plotted as the 
reflection radius versus dimer kinetic energy, Ek, in dimensionless units, kbR{Ek) versus 
Ek/\Eb\- The results are plotted at several different mass ratios, m2jvri\. As expected 
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the reflection radius decreases monotonically with increasing energy. One interesting feature 
is the dependence on the mass ratio. For larger mijm\ the reflection radius is larger but 
decreases quickly with increasing energy. This indicates that a bound state with constituent 
particle masses that are very different will behave like a large but soft deformable ball upon 
contact with the wall boundary. 

When m2 /mi = 1 in one dimension the problem is integrable and exactly solvable via 
the Bethe Ansatz. From the Bethe Ansatz we get 

pcot 5(p) = — 2kb, (66) 

k b R{E k ) = l]^}^ 1 (67) 
As shown in Fig. [6j the lattice results agree with the solution given by the Bethe Ansatz. 
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FIG. 6: One-dimensional lattice results for the reflection radius versus dimer kinetic energy in 
the shallow-binding limit. Also shown are Bethe Ansatz results for 7712/ m\ = 1 and first- and 
second-order results for the expansion of kbclr described in the text. 

Fig. [6] also shows the first- and second-order analytic results for the expansion of kb^r 
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m2 /mi 


K B a R 




k%Pr 


1 


0.4999(2) 


0.005(7) 


0.002(3) 


2 


0.6065(2) 


-0.074(2) 


-0.006(2) 


4 


0.8747(2) 


0.115(2) 


0.006(2) 


8 


1.2149(2) 


0.460(2) 


0.008(2) 



TABLE I: Coefficients of the effective range expansion for one-dimensional dimer-wall 

scattering. 

derived from our general effective theory. 

Table [I] presents the coefficients of the effective range expansion for one-dimensional 
dimer-wall reflection. The error estimates are from the least-squares fitting used in the 



lattice extrapolation and the effective range expansion. From Eq. (66), we see that the 
Bethe Ansatz gives kbo>r — 1/2 for m 2 /m 1 = 1, with all other coefficients equal to zero. 
This completely agrees with the lattice results. These results are universal and could be 
verified using other theoretical methods or experiments such as cold atomic dimers confined 
to a one-dimensional optical lattice with sharp boundaries. 



2. Numerical Results in Two Dimensions 

Fig. [7] shows the reflection radius versus dimer kinetic energy for the two-dimensional 
system in the shallow-binding limit. Again the results are presented in dimensionless com- 
binations, k,bR{Ek) versus Ek/\Eb\- The reflection radius is somewhat smaller than in 
the one-dimensional case. This is reasonable considering the difference between compression 
of a one-dimensional object versus compression of a two-dimensional object along just one 
dimension. In the first case, we expect the object to resist compression more since it does 
not have another direction to expand along as it is compressed. As in the one-dimensional 
case, we see the same dependence on the mass ratio. At large m^jvnx the reflection radius is 
large at small energies while becoming substantially smaller with increasing energy. Fig. [7] 
also shows the first- and second-order results for the asymptotic expansion of kb^r derived 
from the general theory for shallow two-body bound state reflection. 

The coefficients of the effective range expansion for the two-dimensional system are given 
in Table [TT1 
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FIG. 7: Two-dimensional lattice results for the reflection radius vs dimer kinetic energy in the 
shallow binding-energy limit. Also shown are first- and second- order results for the expansion of 
hbclr described later in the text. 





K B an 


KBTR 


k 3 b Pr 


1 


0.407(5) 


-0.10(28) 


0.00(11) 


2 


0.494(5) 


-0.18(19) 


-0.01(7) 


4 


0.718(6) 


-0.01(10) 


-0.00(3) 


8 


1.009(9) 


0.29(6) 


-0.001(17) 



TABLE II: Coefficients of the effective range expansion for two-dimensional dimer-wall 

scattering 

3. Numerical Results in Three Dimensions 



Fig. [8] shows the reflection radius versus dimer kinetic energy for the three-dimensional 
system in the shallow-binding limit. Again the results are presented in dimensionless combi- 
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FIG. 8: Three-dimensional lattice results for the reflection radius versus dimer kinetic energy in 
the shallow-binding limit. Also shown are first- and second-order results for kbclr described in 
the text. 

nations, kbR{Ek) versus Eg/ \Eb\- The results are plotted at several different mass ratios, 
m 2 /m\. Again the reflection radius is smaller than in the two-dimensional case and nearly a 
factor of two smaller than in the one-dimensional case. Just as in the two-dimensional case, 
this conforms to our expectation that the three-dimensional object will resist compression 
even less than the two dimensional object, since it has two unconstrained dimensions it 
can expand along rather than just one. As in the one- and two-dimensional cases, we once 
again find a direct relationship between mass ratio and the derivative of the reflection ra- 
dius with respect to the dimer kinetic energy. Fig. [8] also shows the first- and second-order 
results for the asymptotic expansion of kbcir derived from the general theory for shallow 
two-body bound state reflection. The coefficients of the effective range expansion for the 



three-dimensional system are given in Table III 



We now address what happens in the limit m^jtrix — > 00. Consider the limit rri2 — > 00 
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m2/m\ 








1 


0.279(8) 


— 0.38(9) 


-0.02(2) 


2 


0.342(8) 


— 0.44(6) 


-0.02(2) 


4 


0.506(8) 


-0.27(4) 


-0.02(1) 


8 


0.731(8) 


0.03(2) 


-0.005(5) 



TABLE III: Coefficients of the effective range expansion for three-dimensional dimer-wall 

scattering. 



L 


p(L) 


R[p(L)] 


11 


8 fin 


81(9) MeV 


2.1(4) fm 


9. 


3 fm 


97(10) MeV 


1.6(3) fm 


7. 


3 fm 


118(10) MeV 


1.3(2) fm 



TABLE IV: Momenta and reflection radii for an alpha particle confined to a cube of length 

L. 

with m\ held fixed. In this limit the effective potential converges to a non-vanishing finite- 
valued function, while the mass of the dimer grows with m 2 . One can check this explicitly 



for the expressions in Eq. (46) and Eq. (62). Given the exponential tail of the effective 



potential, the reflection radius near threshold has a logarithmic dependence on m^jmx. We 
see this behavior in each of the plots in Fig. [6j [7l and [8] 



4- Alpha Particle 



Using lattice effective field theory, we have calculated the energy of an alpha particle in 
a cubic box at leading order in chiral effective field theory at lattice spacing a = 1.97 fm 
for cubic boxes of lengths L = 11.8 fm, 9.9 fm, and 7.9 fm. We use the same lattice action, 
algorithms, and codes as in Ref. [TB]. A review of lattice effective field theory can be found 
in Ref. [23]. The vanishing Dirichlet boundaries are implemented as an essentially infinite 
potential energy at the wall boundaries. We note that the ultraviolet divergences are inde- 
pendent of long-distance boundary conditions, and so no new renormalization counterterms 
are needed. In Fig. [9] we plot the energy expectation value of the alpha particle as a function 
of the Euclidean time projection. The exponential curves shown are best fits to a functional 
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FIG. 9: Energy of an alpha particle confined to cubic box of lengths 11.8, 9.9, and 7.9 fm. 



form 

E{t) = E + ce 

These capture the asymptotic behavior at large t. 



-AEt 



(68) 



Using the asymptotic form in Eq. (68 ) we extract the ground state energy E for the given 



boundary conditions. Using Eqs. ^ and ([3]), we can calculate P(L) and R{P(L)\. The re- 
sults are shown in Table HVl The error bars in Table HVl are one standard deviation estimates 
which include both Monte Carlo statistical errors and uncertainties due to extrapolation at 
large Euclidean time. At leading order we find the root-mean-square (RMS) matter radius 



for point-like nucleons to be 1.53(4) fm. Comparing the reflection radii in Table IV to the 



RMS matter radius of the alpha particle, we find that at low momenta the reflection radius 
of the alpha particle is larger than the RMS matter radius. However for larger momenta we 
see a substantial decrease in the alpha particle reflection radius. This suggests that the alpha 
particle is quite compressible under confinement pressure. This appears consistent with the 
observation that alpha clusters are compressed in size when confined within nuclei. Much 
more numerical work is planned to study alpha particles and other nuclei under confinement 
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pressure. 



V. SUMMARY AND FUTURE WORK 

In this paper we derived the phase shift due to the scattering of a two-particle bound 
state on a hard wall in the adiabatic limit and up to second order in an expansion of the 
effective Hamiltonian. We have presented the effective Hamiltonian as an expansion in the 
parameters e~ KBr+ and e~ KBT ~ . We plotted the effective potential in one, two, and three 
dimensions, and discussed the convergence of the expansion. For the equal mass case in one 
dimension we have presented an exact analytic solution using the Bethe Ansatz method. 
Additionaly, we presented numerical and effective field theory calculations of the reflection 
radius in one, two, and three spatial dimensions. We then discussed the consistency among 
the analytic, numeric, and effective field theory calculations. We saw that, as expected, 
the bound state becomes more readily compressible as the number of spatial dimensions 
increases. Furthermore, we noted an interesting mass ratio dependence of the reflection 
radius, noting that for larger mass ratios the reflection radius is larger but decreases faster 
as the kinetic energy of the bound state increases. Using lattice effective field theory we 
calculated the alpha particle energy in a cubic box for L = 11.8 fm, 9.9 fm, and 7.9 fm. We 
then calculated the alpha particle reflection radius and in our analysis of the alpha particle 
compressibility we found that the alpha particle actually appears to compress rather easily. 

We have discussed many aspects of the elastic scattering of quantum bound states from 
a hard surface, paying particular attention to universal behavior which is common to many 
different systems. There appear to be many applications, ranging from experimental predic- 
tions for quantum dots and wells to numerical calculations of nuclear structure and elastic 
deformation. The theoretical techniques employed in this analysis may prove useful in de- 
scribing the effective field theory of other inhomogeneous systems. One immediate extension 
of our three-dimensional analysis is to calculate the deuteron reflection phase shift from ab 
initio lattice chiral effective field theory to verify the universal effective range expansion 



coefficients in Table III, as well as to measure the spin dependence of the reflection phase 
shift resulting from the D-wave component of the deuteron wavefunction. Some other in- 
teresting extensions include investigating the inelastic threshold, and varying the boundary 
conditions used. 
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Another immediate extension is to verify the logarithmic m 2 /mi threshold dependence 
in various physical systems. We expect this effect to be prominent for halo nuclei with a 
heavy core and one satellite nucleon as well as in quantum dots and wells for semiconductors 
with a large ratio between hole and electron effective masses. It could also be reproduced 
experimentally with heterogeneous cold atomic dimers consisting of one heavy and one light 
alkali atom. 

Financial support from the U.S. Department of Energy and NC State G A ANN Fellowship 
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Forschungszentrum Julich. 



[1] T. Takagahara and K. Takeda, Phys. Rev. B 46, 15578 (1992). 

[2] J. Y. Marzin, J. M. Gerard, A. Izrael, D. Barrier, and G. Bastard, Phys. Rev. Lett. 73, 716 
(1994). 

[3] S. Pfalz, R. Winkler, T. Nowitzki, D. Renter, A. D. Wieck, D. Hagele, and M. Oestreich, 

Phys. Rev. B 71, 165305 (2005). 
[4] L. He, G. Bester, and A. Zunger, Phys. Rev. Lett. 94, 016801 (2005). 

[5] S. Montangero, D. Frustaglia, T. Calarco, and R. Razio, Europhys. Lett. 88, 30006 (2009). 
[6] M. Liischer, Commun. Math. Phys. 104, 177 (1986). 
[7] M. Liischer, Commun. Math. Phys. 105, 153 (1986). 

[8] S. R. Beane, P. F. Bedaque, K. Orginos, and M. J. Savage (NPLQCD), Phys. Rev. D73, 

054503 (2006), hep-lat/0506013. 
[9] S. R. Beane, P. F. Bedaque, K. Orginos, and M. J. Savage, Phys. Rev. Lett. 97, 012001 (2006), 

hep-lat/0602010. 

[10] V. Bernard, M. Lage, U.-G. Meifiner, and A. Rusetsky, JHEP 08, 024 (2008), arXiv:0806.4495 
[hep-lat] . 

[11] E. Epelbaum, H. Krebs, D. Lee, and U.-G. Meifiner, Phys. Rev. Lett. 104, 142501 (2010), 

arXiv:0912.4195 [nucl-th]. 
[12] E. Epelbaum, H. Krebs, D. Lee, and U.-G. Meifiner, Eur. Phys. J. A45, 335 (2010), 

arXiv: 1003.5697 [nucl-th]. 
[13] A. Tohsaki, H. Horiuchi, P. Schuck, and G. Ropke, Phys. Rev. Lett. 87, 192501 (2001), nucl- 



28 



th/0110014. 

[14] M. Chernykh, H. Feldmeier, T. Neff, P. von Neumann-Cosel, and A. Richter, Phys. Rev. Lett. 

98, 032501 (2007). 
[15] Y. Suzuki et al., Phys. Lett. B659, 160 (2008), nucl-th/0703001. 

[16] E. Epelbaum, H. Krebs, D. Lee, and U.-G. MeiBner, Phys. Rev. Lett. 106, 192501 (2011), 

arXiv:1101.2547 [nucl-th]. 
[17] D. Lee and M. Pine, Eur. Phys. J. A47, 41 (2011), 1008.5187. 
[18] S. Kdnig, D. Lee, and H.-W. Hammer, Annals Phys. 327, 1450 (2012), 1109.4577. 
[19] S. Kdnig, D. Lee, and H.-W. Hammer, Phys. Rev. Lett. 107, 112001 (2011), 1103.4468. 
[20] S. Tan (2012), 1201.6014vl. 
[21] J. Macek, J. Phys. B 1, 831 (1968). 
[22] C. D. Lin, Phys. Rev. 257, 1 (1995). 

[23] B. D. Esry, C. D. Lin, and C. H. Greene, Phys. Rev. A 54, 394 (1996). 
[24] D. Lee, Prog. Part. Nucl. Phys. 63, 117 (2009), arXiv:0804.3501 [nucl-th]. 

Appendix A: Derivation of the Diagonal Adiabatic Correction 

We have already defined the diagonal adiabatic correction as 



T(X) 



1 



\1>x), 



(Al) 



2M 



dX 2 



where 



(A2) 



The normalization condition gives 




) fd ( K r_,r_ 




(A3) 



where 
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(t>r-,r+{r) = Vdfi{^r-,r+, \A) ~ ^,o(^r_,r+, \f ~ f-\) ~ yd,o(«r_,r+, \r ~ f+\) 

+ yd,o(«f_,f+, \r - f+ + f_|) + yd,o( K f-,r+,r + f+ - f_) H 



(A4) 



and A?_f + is the function that normalizes the relative coordinate wavefunction, (f)?_f + (r). 
The ci-dimensional Yukawa function is 



{2n) d p 2 + k 2 



(A5) 



and the overlap integral is 



/„ («W) = | rfV ydfi (iJ, \r'\)y d ,o(K H , \r'-r\), (A6) 

for two Yukawa functions whose centers are separated by a spatial distance r. 

We think of </>?_,?+ (0 as a Yukawa function centered at the origin, plus images from a 
single mirror reflection, plus images from two mirror reflections, and so on. We expand 
-*4f_,F+ by including more reflections at each increasing order. This procedure gives 



A 



r- ,f+ (0) 



[fd(K B ,K B ,0)] 1/2, 



and 



(A7) 



. 1 fd (kb, *b, r_) + f d (k b , k b , r+) 

A r-,r+(l) — 9 3/2 + «■ 

2 [/d(KB,KB,0)] 



<9A 



f_,f+(0) 



1 fd(K B ,K B ,r + ) + / d (/Cfl,K B ,r_) - K?_,? + (i)Q^f d (kb,k b ,0) 
~2 [/d(«B,« B ,0)] 3/2 ' 

Since the relative coordinate wavefunction is normalized, 



(A8) 



2^*1 ^2 l« 



2 MAI - 

r_,r+ 



+ 



2M 



1 d d ~ l r\ 



<9X 



(r) 



(A9) 



We compute this term as an expansion in powers of e Ksr + and e KBT - . At leading order 
we find 
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1 d 2 

(lpX\ -T^TT \^x) 



2M XT ^' dX d 
M 

2mlf d (K B ,K B ,0) 



(i) 



Vdjo («b, r+) + 



M 



2m5/«j «b, 0) 



2/<i,o («s,r- 



(A10) 



At second order 



1 



O 2 



(2) 



2M xr "' dX 2 

T^(X)+T^(X)+T^(X) 



(By 



+ tSax) + t% 2) {x) + rg> (X) + r« (x) + r<g(x), 



(All) 



where the first term in Eq.(All) is 



T§{X) 



fd («B, «B, 0) / 9v4 r -_ ir - + (i) 



2M 



<9X 



+ 



(A12) 



and the second term 



Tll](X) 



M 



r + + r_ 



-yd,-\(K>B,r+ + r. 



(A13) 



l~l 2 fd(KB,KB,0) ' ' /i/d(«B,«S,0)' 

contains terms due to two derivatives with respect to the spatial coordinate without any 
contribution from the first-order corrections A?_f + m and Kf-.f+fl)- 



Af/t 



r_,f+(l) 



5 



2mf/ rf (k, «, 0) 9k 



+ 



r_,r+(l) 



<9 



2mlf d (K,K,0) 8k 



+ 



.4 



F_,r+(1) 



[/d(«B,«B, 



,1/2 



M M . 
—2^,0 (kb, r + J H 2^>o («b, r_) 



n? 



contains terms proportional to «?_,?+(!) and Aj? _,r + (i)- 
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(Dl) 
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r_,r + (l) 



mi/d (k b , ^b, 0) <9X 



(«B,r + ) + yy d)2 (KB,r H 
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, (A15) 
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tIU x ) 
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m 2 f d {k b , k b , 0) dX 



V T 3 
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Ok, 



r-,r+(l) 



Amxfd (k b , k b , 0) dX 



r + y d ,2 (k b , r+) + 3y d)3 (k b> r+)] 



and 



Kb 



8k. 



f_,f+(i) 



[r-yd,a («, r_) + 3y d>3 (re, r_)] 



, (A16) 
(A17) 

(A18) 



contain terms due to a derivative with respect to k as well as a spatial derivative. Finally, 



dX 



dy, 



3- — («b, 0) + 3k 2 (k b ,0) - k 3 (kb,Q) 



d 3 y d ,o 



8k% 



24MK B f d (K B ,K B ,0) [ ~ 8n B 
contains terms due to two derivatives with respect to k. 



8k% 



(A19) 



Appendix B: Diagonal Adiabatic Correction in One Dimension 



Using Eq. (A5) and Eq. (A6) we compute the adiabatic correction in one spatial dimen- 



sion. The first term in Eq. (All) 



e -2K B r +K 2 B + m ^ J e K B (r++r_) TOl (3 + Kgr _) + m2 (3 _ KbT+ 

2m\m\ 



(Bl) 



and the second term in Eq.(All), due to two spatial derivatives but not including contribu- 



tions due to Ap_ t ? + n) and K?_ ) ? + m, is 



2 e ~ K B(r++r-) _|_ m ^ ^_ m 2 _ m 2 _|_ mim<2 ^ B r + + K B T — — 2]) 



2 2 

m\m2 



(B2) 



The term proportional to Kf_ ,♦=+(!) an d ^4r_,r+(i) is 
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= 5 — k 2 b (mi + m 2 ) (2k b t_ - 1) H 5 — 4 ( m i + m 2) (2K B r + - 1) 

m 2 mf 
e - Ks (r++r_) / 2 + m 2) 

+ (mi + m 2 ) (2« B r + + 2« B r_ - 1) . (B3) 



mfm 2 

The four terms due to one derivative with respect to k and one spatial derivative are 



T$ X) {X) = (B4) 
2m 2 m| 12 + 1 j ' 



^2)W = 

k% (mi + m 2 ) (9 - AK B r + + n B r 2 _) 
2mim 2 

T^JX) = " l(mi + m 2 2)(3 + ^ r+) (m 2 e-— + +mie--^-)) , 
y A > mfm 2 v ' 



(m 2 e- Kfl(r++r -) + mie- 2Ksr -) , (B5) 



(B6) 



and 



(2) (x) = 4K+m 2 )(3 + «Br-) + mie - 2 « s ,-) . (B7 ) 

1 j mim 2 

The term containing contributions due to two derivatives with respect to k is 



p ~2K B r + -2/t s r_\ i n -ft B (r + +r_ ) 2 

7g' ) (X) = 5M4(^ r + ^ r ) + 1 ° e ' W 



Appendix C: Diagonal Adiabatic Correction in Two Dimensions 



Using Eq. (A5) and Eq. (A6) we compute the adiabatic correction in two spatial dimen- 
sions. We find that in two dimensions the first order correction to the binding momentum 
is 

«r_,r+(i) = -k b {K (n B r + ) + Kq(k b T-)\ , (CI) 
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the second order correction to the binding momentum is 



«F_,f+(2) = -y {^o(^+) 2 + K (K, B r_) 2 + AK [k b (r+ + r_)]} 

- K B K (K B r + ) [n B r + Ki(K B r + ) + K B r— K\ {k b T— ) - if («B7--)] 

- K 2 B Ko{K B rJ) \r + K 1 (K B r + ) + « B r_lifi(KBr_)] , 



(C2) 



and the generalized Yukawa function in two dimensions is 

y2,o( K > r ) = 7T K o(Kr). 
Zrr 



(C3) 



The first term in Eq. (All) 



%m\m\ 



{K B mir_Ks(K B rJ) + K B m 2 r + K G (K B r + ) - 67^1^ 



-Qm 2 K l (K B r + ) + K B m x r_K 2 {n B rJ) + K B m 2 r + K 2 (K B r + )} 2 , 



(C4) 



and the second term due to two spatial derivatives but not including contributions due to 

^4f_,f+(i) and /tr_,r + (i), is 



2Mk| 

2 2" 
mfm 2 



{-(mi + m 2 ) 2 K [KB(^+ + r_)] 



+K B miin 2 (r + + r_)i^i[K B (r + + r_)]} . 



(C5) 



The term proportional to /v_,r + (i) and A*? _,f + (i) is 



mi 



[K Q {Kr + )\ 



K*M[K {Kr + ) + K Q {Kr.)] d 
+ ml dK 



+ 2/tv / vr^_/ + (i) 



M 
2irml 



K (K B r^ 



M 
2irm 2 



K (K B r. 



(C6) 



where 
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Ar-,r+{i) = V^kb [-2K (K B r-) - 2K {n B r + ) + K B r-Ki(K B r-) + K B r + K 1 (K B r + ) 



-1/2 



(C7) 



The four terms due to one derivative with respect to k and one spatial derivative are 



Mk 2 r , . . I - 2 2 / .2.1 



4m|m 2 



[mi^i(K B r-) +m 2 ifi(K B r + )] < Sk b t + G 1 



K B r + 1 





2 ' 2 


_I I _i 




2' 2' 



+3tt [/t jB r + K (/« J B^+) + 2K 1 (K B r + )\ Li(/c B r+) + K B r + K 1 (n B r + ) [An B r + + 37rL 2 (/t jB r + )]} , 

(C8) 

where L a is the modified Stuve function and is the generalized Meijer G-function. 



\m x K x {K B r_) +m 2 K 1 (K B r + )] t -?>K B r 2 _G 2 l 'l 



K B r^ 1 





2 ' 2 


_I I _i 




2' 2' 



+37r [KBr-KoC^Br-) + 2Xi(/« B r_)] Li(/c B r_) + K B r_Ki(K B r^) [Ak b t^ + 37rL 2 (/t J Br_)]} . 

(C9) 



B <9/V_,r+(l) 



K 2/2,2 («b, r+) + 3y 2 ,3 («b, r+)] , 



(CIO) 



where the spatial derivative of the first order correction to the binding momentum is 

oX mi m 2 

The generalized Yukawa function 2/2,2 r ) is 



(Cll) 



K (Kr)L^i(Kr) — Ki(nr)L 



27TK 



tfi(«r), (C12) 
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the generalized Yukawa function y 2 , 3 (ft, r) is 



2/2,3 



nr 



— K (K,r)L_i(K,r) — Ki(K,r)L (nr) 



+ 



1 



247TK 3 



Ki(nr) [—AK 2 r 2 + 67tLi(kx) + 37r/trL 2 (/tr)] 



J_ q2,1 

8ttk 1,3 







kx 1 





T' 2 


_1 I _1 




2' 2' 



and 



TTK, 



8k 



r_,f+(l) 



[r_y 2 ,2(K,r_) + 3y 2 , 3 (ft, r_)] . 



( D4 ) V "' m 2 «9X 
The term containing contributions due to two derivatives with respect to k is 



^(2) / v \ 8k b ( m i + m 2) [miii'i(K B r_) + m 2 iTi(/t B r + )]^ 



Appendix D: Diagonal Adiabatic Correction in Three Dimensions 



(C13) 



(C14) 



(C15) 



For the effective Hamiltonian in three dimensions /Cf_,r+(i)j K r-,r+(2), T^(X), T^(X), 
T^(X), T^ 1D4 ^(X), and T^(X) are defined as follows. The first order correction to the 
binding momentum is 



K 



r-,r+(l) — 



J/3,o(«b,?+) +y 3fi (K B ,r-) 



8k, 



; 1/3,0 («b, 0) 



(Dl) 



the second order correction to the binding momentum is 



K 



r-,r+(2) 



d 
8k B 



2/3,0 («B, 0) 



2/3,0 («B, 7+)+ 2/3,0 (ftB, 7"-) g 
^2/3,0 («B, 0) 0«B 



[2/3,0 («B,r+) +J/3.0 («B,r_)] 



-2y3,o(«B,r+ + r_) - 



[2/3,0 («B,r+) +2/3,0 («B,r_)] 2 <9 2 



d 2/3,0 («b,0) 



9k s 



94 



S/3,o(«B,0) > , 



(D2) 
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and the generalized Yukawa function in three spatial dimensions is 



47rr 



The first term in Eq. (All) is 



(D3) 



(A) 



h ( k b, Kb, 0) / dA ? _f + (i) 



2M V OX 
where the overlap integral in three spatial dimensions is 



(D4) 



is («, «, r) 



87TK 



(D5) 



and the spatial derivative of the first order correction to the wavefunction normalization is 

C^4r_,r+(l) 



dX 

2M a 



1 if ^h( K B,KB,r + ) + ^^-f 3 {K B ,K B ,r- 



2 [fs(KB,KB,0)f 2 

The next term, with contributions due solely to two spatial derivatives, is 



(D6) 



y 3 ,o{K B ,r + + r„) + ——± ^— 2/3 _!(/«£,, r+ + r_), (D7) 



A* 2 /3 («B, «s, 0) ' ' /i/ 3 («B, Kb, 0) 

where the generalized Yukawa function, y 3) _i(AC, r), is 



J/3,-i(«,r) 



5 



The term containing contributions due to Kf_,f+(i) an d ^4r_,r+(i) i 



is 



(D8) 



(C) 



Mk 



f— ,r+ (1) 



5 



2m\fz (k, k, 0) 9k 



+ 



Mk 



r_,r+(l) 



<9 



1m\f 3 («, «, 0) 9k 



[2/3,0 («, ?+)] 



[2/3,0 («,r_)] 



K=Kb 



+ 



A 



r_,r+(l) 



[/3(«B,KB,' 



,1/2 



M 



M 



r 2/3,0 H 2^/3,0 (,«B,r 



(D9) 
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The four terms that result from taking one derivative with respect to k and one spatial 
derivative are 



k b 



8k? 



F-,P+(1) 



T T' 3 

y?/3,l( K B,r + ) + yt/3,2(KB,r + ) + -t/ 3 ,3(KB,r + ) , (D10) 



m 1 f 3 (K B ,K B ,0) 8X 
where the spatial derivative of the first order correction to K?_f + is 



8k? 



?-,?+{!) 
8X 



2M 



J_ 9 

mi 



9^2/3,0 («B,r+) + ^^J/ 3l o(«B,r_) 



9^2/3,0 («B,0) 



and the generalized Yukawa functions, 2/3,1 (k, r), 2/ 3 , 2 (/t, r), and yz, 3 {K, r), are 



(Dll) 



and 



J/3,1 = -— Ei(-«r), 
M«,r) = ^ + £ F Ei(-«r), 



(D12) 
(D13) 

(D14) 



Here Ei is the exponential integral function. 



tIU x ) 



k b 



8k 



r_,r+(l) 



m 2 h («b,«b,0) <9X 



r r_ 3 

yS/3,1 («B, r_) + yl/3,2 («b, r_) + "2/3,3 («b, r_) 



<9k. 



f_ ,r+ (1) 



and 



4mi/ 3 (k b , «b, 0) <9X 

r (2) m K B ^r_,r + (l) 

( S4 ) 1 J 4m 2 / 3 (« B ,« B) 0) 8X 



[r+2/3,2 (kb, r+) + 3?/ 3 ,3 («b, , 



[r- 2/3,2 («,r_) + 32/3,3 («,r_)]. 



, (D15) 
(D16) 

(D17) 



The term containing contributions due to two derivatives with respect to k is 



( ^r_,r + (l)\' 
I 9X J 



(E) 



24MK 3 B f 3 {K B ,K B ,0) 



0^2/3,0, n , . q <9 2 2/ 3 ,o 
—3 — — (kb, Uj + 3/t- 



8 k 1 



8k\ 



k b , 0) - k 



2 d 3 2/3,o 
8k% 



K B ,0) 



(D18) 
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Appendix E: Summary of Continuum Extrapolations 



In one dimension for mass ratio mijm\ = 1 we consider L z \c\ = 8, 10, 20, and for L z \c\ = 
10 we consider the ground state as well as the first and second excited states. Table V shows 
a summary of the continuum extrapolations in one dimension for mass ratio mijmx = 1. 



L z \c\ 


L z 


Rkb 


Er/\Eb\ 


20 


200, 300, 400 


0.498 


0.006832 


10 


100,200,300,400 


0.4950 


0.030 


8 


80,160,240 


0.4981 


0.050122 


10 


100, 200, 300 


0.4812 


0.12083 


10 


100, 200, 300 


0.461 


0.2694 



TABLE V: Summary of continuum extrapolations for mijmx = 1 in one dimension 



Table VI gives the values of L z \c\ and L z considered as well as the continuum limit 



extrapolations in one dimension for mass ratio mijmx = 2. 



L z \c\ 


L z 


Rub 


Er/\Eb\ 


20 


200, 300, 400 


0.6041 


0.0062 


10 


100, 200, 300 


0.5957 


0.028228 


8 


80, 160, 240 


0.5885 


0.0471 


10 


100, 200, 300 


0.5676 


0.111164 


10 


100, 200, 300 


0.5311 


0.247 



TABLE VI: Summary of continuum extrapolations for m 2 /mi = 2 in one dimension 



In one dimension for mass ratio mijmx = 4 the continuum extrapolations are given by 
Table MR 
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L z\c\ 


L z 


Rkb 


Er/\Eb\ 


20 


200, 300, 400 


0.8694 


0.00473 


10 


100, 200, 300 


0.8504 


0.02292 


8 


80,160,240 


0.8343 


0.039 


10 


100, 200, 300 


0.7915 


0.0891 


10 


100, 200, 300 


0.7214 


0.19409 



TABLE VII: Summary of continuum extrapolations for mj/m] = 4 in one dimension 

In one dimension for mass ratio 1712/mi = 8 we get the continuum extrapolations as 
shown in Table IVIII1 



L z\c\ 


L z 


Rkb 


Er/\Eb\ 


20 


200, 300, 400 


1.206 


0.003151 


10 


100, 200, 300 


1.1742 


0.01665 


8 


80, 160, 240 


1.1455 


0.029907 


10 


100, 200, 300 


1.083 


0.0635 


10 


100, 200, 300 


0.982 


0.13589 



TABLE VIII: Summary of continuum extrapolations for m 2 /m 1 = 8 in one dimension 

In two dimensions for mass ratios 7712/ mi = 1,2,4,8 we consider the ground and first 
excited state for each L z kb = 6.8, 10.2. As noted previously L z is the length of the lattice 



in the direction perpendicular to the confining wall. Table IX shows a summary of the 
continuum extrapolations in two dimensions for mass ratio m^jvnx = 1. 



L Z K B 


L z 


Rkb 


Ek/\Eb\ 


10.2 


60,90,120,150 


0.485146 


0.004742 


6.8 


40,60,80,100 


0.456299 


0.011126 


10.2 


60,90,120,150 


0.482144 


0.019341 


6.8 


40,60,80,100 


0.45136 


0.044811 



TABLE IX: Summary of continuum extrapolations for m 2 /mi = 1 in two dimensions 

Table [X] gives the values of L z k,b and L z considered as well as the continuum limit 
extrapolations in two dimensions for mass ratio m^jvrLx = 2. 
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L Z K B 


L z 


Rkb 


Er/\Eb\ 


10.2 


60,90,120,150 


0.572157 


0.004342 


6.8 


40,60,80,100 


0.541738 


0.010386 


10.2 


60,90,120,150 


0.567111 


0.017369 


6.8 


40,60,80,100 


0.532998 


0.040117 



TABLE X: Summary of continuum extrapolations for mijm\ = 2 in two dimensions 

In two dimensions for mass ratio m^jvax = 4 the continuum extrapolations are given by 
Table E0 



L Z K B 


L z 


Rkb 


Er/\Eb\ 


10.2 


60,90,120,150 


0.791197 


0.003473 


6.8 


40,60,80,100 


0.759959 


0.007583 


10.2 


60,90,120,150 


0.783424 


0.013027 


6.8 


40,60,80,100 


0.742787 


0.031096 



TABLE XI: Summary of continuum extrapolations for m 2 /wi = 4 in two dimensions 

In two dimensions for mass ratio m 2 /m 1 = 8 we get the continuum extrapolations as 
shown in Table IXII1 



L z kb 


L z 


Rkb 


E K /\E B \ 


10.2 


60,90,120,150 


1.81456 


0.09496 


6.8 


40,60,80, 100 


1.1018 


0.17597 


10.2 


60,90,120,150 


1.21421 


0.20091 


6.8 


40,60,80,100 


0.79017 


0.43994 



TABLE XII: Summary of continuum extrapolations for 7772/7771 = 8 in two dimensions 
In three dimensions we consider the ground and first excited states of L z kb = 4.9,6. 



Table XIII shows a summary of the continuum extrapolations in three dimensions for mass 
ratio 777,2/777.1 = 1 



Table XIV gives the values of L z kb and L z considered as well as the continuum limit 
extrapolations in three dimensions for mass ratio m^jvrii = 2. 



41 



L Z K B 


L z 


Rkb 


Ek/\Eb\ 


6 


24,36,48 


0.27263 


0.08294 


4.9 


20,30,40 


0.26891 


0.12394 


6 


24,36,48 


0.25605 


0.32787 


4.94 


20,30,40 


0.24773 


0.48669 



TABLE XIII: Summary of continuum extrapolations for m^jmi = 1 in three dimensions 



L Z K B 


L z 


Rub 


ErI \Eb\ 


6 


24,36,48 


0.32937 


0.0769 


4.9 


20,30,40 


0.32308 


0.11573 


6 


24,36,48 


0.30134 


0.30132 


4.9 


20,30,40 


0.28801 


0.44853 



TABLE XIV: Summary of continuum extrapolations for 777.2/ 'mi = 2 in three dimensions 
In three dimensions for mass ratio mijmi = 4 the continuum extrapolations are given 



by Table XV 



L Z K B 


L z 


Rub 


ErI \Eb\ 


6 


24,36,48 


0.47906 


0.06215 


4.9 


20,30,40 


0.46579 


0.09545 


6 


24,36,48 


0.42392 


0.23812 


4.9 


20,30,40 


0.39833 


0.35779 



TABLE XV: Summary of continuum extrapolations for m^jmi = 4 in three dimensions 

In three dimensions for mass ratio m^jmx = 8 we get the continuum extrapolations as 
shown in Table IX VII 
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L Z K B 


L z 


Rk>b 


Er/\Eb\ 


6 


24,36,48 


0.68252 


0.0454 


4.9 


20,30,40 


0.65848 


0.07191 


6 


24,36,48 


0.59503 


0.16868 


4.9 


20,30,40 


0.55407 


0.25769 



TABLE XVI: Summary of continuum extrapolations for 1712/1711 = 8 in three dimensions 
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